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Abstract 

We consider the high energy behaviour of the amphtudes for pair production of charged leptons, 
quarks, Higgs bosons, sleptons, squarks and charginos at lepton colUders. We give the general 
expressions of the leading quadratic and subleading linear logarithms that appear at the one loop 
level, and derive the corresponding resummed expansions to subleading logarithmic order accuracy. 
Under the assumption of a relatively light SUSY scenario and choosing the MSSM as a specific 
model, we compare the predictions of the one-loop and of the resummed expansions at variable 
energy. We show that the two predictions are very close in the one TeV regime, but drastically 
differ in the few (2, 3) TeV range. 
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I. INTRODUCTION 



It is by now well-known that the electroweak radiative corrections to Standard Model 
pair production processes increase strongly with the center of mass energy -^i at the one- 
loop level. This is due to the presence of large double and single logarithms ~ - log^ -rp-, 

■ ■ ■ ■ ■ ■ ■ ■ ■ ■ w 



. In the TeV range such terms reach the several percent level, which 
should be easily observable (and measurable) at future lepton colliders 0, Q]- Actually, for 
energies beyond the few TeV range, the numerical size of their effect begins to be too large 
(beyond the relative ten percent level), and the validity of the simple one- loop approxi- 
mation must be seriously questioned. This has led a number of theoretical groups to propose 



resummation prescriptions 



Without entering the details of the different calcu- 



lations, we shall accept the conclusion that, for massive final pairs, an agreement seems to 
exist that a full resummation can be given, but only to subleading logarithmic level. 
In the case of supersymmetric extensions of the SM, similar studies have been recently 
performed for the processes of production of fermion pairs at the one-loop level and 
of scalar pairs |131] in the MSSM. In ref. both the one-loop approximation and the 
resummed expansions have been computed. Under the qualitative assumption of a rela- 
tively light SUSY scenario, with all the relevant masses of the process " adequately" smaller 
than e.g. one TeV, it has been shown that the two approximations are "essentially" (i.e. 
at the expected relative one percent level of experimental accuracy) identical in the ~ 1 



TeV region (final possible reach of the proposed LC p), but differ drastically in the few 
(2,3) TeV range (aim of the proposed CLIC 7]). Thus, for the spin zero production case, 
to subleading logarithmic accuracy, a simple one-loop approximation would seem definitely 
inadequate at CLIC energies, but valid in the LC regime. 

The aim of this paper is that of performing a general investigation, analogous to that 
already carried through for scalar production, to include also the processes of charged spin 
one-half (fermions, charginos) pair production. For all these cases we shall give both the 
one-loop level and the resummed logarithmic expansions to subleading logarithmic accu- 
racy. Under the assumption of a relatively light SUSY scenario, we shall then compare the 
two approximations at variable cm. energies for a number of experimental observables. 
We shall show that the same conclusions obtained for scalar production are apparently 
valid for the extended case, making the validity to subleading logarithmic order of a one- 
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loop approximation for spin zero and one-half charged pair production in the MSSM at 
future colliders to only depend on the chosen cm. energy, and not on the specific process. 
Our analysis has been performed and our conclusions have been drawn in the particular 
case of the MSSM, but they could be easily adapted or modified to treat different supersym- 
metric models, although we believe that this is beyond the purposes of this first, necessarily 
preliminary work. 

Technically speaking, the plan of this paper is the following: Section 2 contains the 
one-loop expansion of all the considered processes. The resummed expansions are shown 
in Section 3, and a comparison between the two approximations is given in Section 4. A 
final discussion and some possible conclusions are given in Section 5. The definitions of the 
observables for final fermionic or scalar pairs are given in Appendix A. Appendices B and 
C contain a summary of several long (but necessary) analytic formulae for the asymptotic 
expansions at the one loop level and with resummation. 



II. ASYMPTOTIC LOGARITHMIC EXPANSION AT ONE-LOOP 



A. Generalities 



As illustrated and discussed in several previous papers llJ| , at the one- loop level, 

the logarithmic terms appearing in e^e'^ annihilation processes can be separated into three 
categories, Renormalization Group (RG) terms. Universal terms and Non Universal (angular 
dependent) terms, so that the invariant scattering amplitude for the pair production process 

e+e- hh (2.1) 
(a representing the electron chirality) can be written in the following form 

^1 loop ^RG _j_ ^univ _j_ ^non univ ^2 2) 

The RG contribution represents the linear logarithms Q| generated by the "running" of 
the gauge coupling constants, that are known and calculable in a straightforward way. It is 
obtained by introducing in ^4^°™ the running SU{2) x U{1) couplings g^g' according to the 
asymptotic MSSM /3o, /3q functions: 

/^-^^-f-f (2.3) 



As) = —j^ — , 9%) = yyj — (2.4) 

where Ca = 2, = 3, n/i = 2 and gsm'^}^, = gf'cosi?w = e. 

At one loop, this quite general procedure gives single logarithms that do not factorize, 
but can be obtained from A^"^"^ as: 

where // is a reference scale defining the numerical values of g,g'. 



The universal contribution can be of quadratic and of linear kind and in a covariant 
gauge is produced by diagrams of vertex (initial and final triangles) and of box type; it 
only depends on the quantum numbers of the initial e~^, e" and final /i, /2 lines and can be 
written as: 

^univ^^Born_ |^ ^in ^ ^fin j (2.6) 

with the correction to the initial e^e^ lines 



{le, + 1) + 1 I 2 log ^ - V ^ J (2.7) 



where My is a common gauge boson mass (we use Mz ^ Mw for simphcity and we have sep- 
arated the infrared part of the photon contribution, keeping here only the ultraviolet one by 
putting My = Mz), 01 is the chirality and /, Y the isospin and hypercharge of the initial elec- 
tron e~ . The expression for will be given separately for each of the examples listed below. 

The uou universal (angular (lependent) contribution can be written 

^non univ ^Born _ q^^^S, (2 g) 

It only consists in residual terms arising from the quadratic logarithms log^ t, log^ u (from 
which the log^ s part has been subtracted and put in the universal contribution) generated 
by box diagrams containing W,Z,^ gauge boson internal lines where t — — f (1 — cos'j?) 
and u — — I (1 -|- cos{}), being the scattering angle. There are only few such diagrams 
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and they have been all explicitely computed. 



To subleading logarithmic accuracy, these three types of contributions have been calcu- 
lated exactly. We shall now summarize, trying to make the review as short as possible but 
reasonably complete and self-contained, the results of our effort for the charged final pairs 
of the MSSM, including those for scalar production and fermion production that were al- 
ready given in other papers. The list of considered final states is now given in the following 
subsections. 



B. Chiral Lepton or Quark Pair f^fp 

The Born amplitude is given in Appendix Bl, using the notations given in Appendix 
A; for more details see previous papers For a final fermion pair with chirality (3, the 
logarithmic part of the one loop amplitude (j2.2|) is obtained by adding eq. (j2.5|) . (j2.6j) . (j2.8p 
with: 



^fin gauge 



+ C 



fin Yukawa 
/3 



(2.9) 



and 



c 



fin gauge 
'/3 



lib {^'"^'" + 1' 



y2 ■ 



2 log 



log^ 



c 



fin Yukawa 

/3 



'1 + 5/3,R rhj rh), 



167r2 



M2 



/3,L 



2M^ 



(2.10) 
(2.11) 



^ang 



9 1 s 
167r2 °^M2 



H^wV;,>7. + 4/Vt)iog^ 



5dj log — - (5„ J log — 

' s s 



(2.12) 



We denote irif = rrit/ sin/5 if / = t andmj = m^/ cos/5 if / = 6; /' denotes the corresponding 
isopartner of /, and (3 (not to be confused with the chirality index) is the mixing angle 
between the vacuum expectation values of the up and down Higgs chiral superfield (in 
standard notation tan/5 = Vu/vd)- 

For particles other than those belonging to the third family of quarks/squarks, the Yukawa 
terms are negligible. 
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C. Slepton or Squark Pair f^fp 

The results for this case are quite similar to those valid for fermion production, see 
Appendix B2 and ref . . The logarithmic part of the one loop contribution for the process 

e+e- (2.13) 

{j3 = L, R) can be written in the same form as for e'^e'^ — > f pfp producing a fermion pair 
with chirality with the same expressions for c™, c^^ s^^s'^^ (Sm Yukawa ^ang^ 

D. Charged Higgs Bosons or charged Goldstones 

In the MSSM the charged Higgs bosons and Goldstones are produced in pairs through 
the same diagrams that appear in the case of sfermion production. At / s accuracy the 
amplitudes for G^G^ production are equivalent to the physical amplitudes for longitudinal 
W^Wj^ states. The previous equations, (j2.1UII2.l2|) which concern the gauge parts give the 
correct result for the process 

ete-~^H+H- (2.14) 
provided that is replaced by with the following quantum numbers 

Q{H-) = -1, I\H~) = -\ (2.15) 



For what concerns the Yukawa part, one has 

fin Yukawa _ .s" / m^ot^ /3 tau^ /? \ s 

These formulae then apply also to 

e+e-— ^G+G- (2.17) 

without any change concerning the Born term and the one loop gauge terms, but the Yukawa 
part has to be simply modified with ml tan^ (3 ml and mf cot^ (3 mf. 

E. Charginos xfxj 

Charginos are mixtures of gaugino(Wino) and Higgsino components selected respectively 
by the mixing matrix elements Zf^, using the notations of ref.^^. The expressions of 
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the Born terms and the one-loop contributions written in the Appendix B3 reflect clearly 
the properties of these two types of components. 

The Born amplitude involves both Wino components {Zfj^Zf^ and Higgsino components 
(Z^^Z^j) produced in s-channel through 7, Z formation, but it only involves Wino com- 
ponents produced through sneutrino exchange in the w-channel (the contribution of the 
Higgsino component vanishes like the electron mass). 

The RG amplitude is computed using Eq. ()2.5|) . We have written separately its Higgsino 
part. Its Wino part is regrouped with the other one loop Wino contributions because of 
remarkable cancellation properties explained below. 

The universal terms c^"^ generated by the xtxj lines contain Higgsino and Wino compo- 
nents that factorize the corresponding Born term components. We have computed them both 
through one loop diagrams in the 't Hooft-Feynman gauge and through the splitting function 



formalism [17| with addition of Parameter Renormalization terms obtaining an agreement 
that we consider, given the not simple structure of the related formulae, encouraging. We 
observed the following properties. The s-channel terms contain universal corrections from 
both Higgsino and Wino components; their Wino part can be identifled with the sum of 
a splitting function contribution and of a Wino Parameter Renormalization term (the 
RG Wino part computed through eq. ()2.5|) . see details in Appendix B3); its Higgsino part 
contains both a "gauge" and a "Yukawa" part. The w-channel terms only contain a Wino 
part also identiflable with the sum of a x^ splitting function contribution and of a Wino 
Parameter Renormalization term. These properties of the Wino and of the Higgsino parts 
are similar to those observed in the cases of W^W~ and H^H^ production, respectively; 
in particular at one loop we get only a DL term (the SL terms cancel) for the universal 
"gauge" parts both in W^W~ and in Wino pair production. 

The non universal (angular dependent) terms arise from residual terms of 77, 7Z', 
ZZ and WW boxes in the s channel and of single 7, Z, and W boxes in the sneutrino 
exchange channel leading to both 1/u and 1/t terms (see Appendix B3). 

A point that must be noticed is the fact that an important difference exists between 
the chargino pair production and the previous considered processes. This is due to the 
fact that, in the chargino case, the Born approximation already contains typical SUSY 
parameters in the mixing coefficients. Being by definition "bare" quantities, a suitable extra 
renormalization is required. This is a well-known fact, already discussed in previous papers 
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(see e.g. |l8|), and the choice of a convenient renormahzation scheme is essential. In our 
case, at the chosen level of logarithmic accuracy, we can neglect this complication since the 
difference between the bare parameters and the physical ones will be in any case a constant 
term of order a, and will not find place in our expansion (it should be, though, properly 
retained in a more complete next-to next to leading order- analysis). This means that in 
our one-loop expansions one can systematically assume, for the values of these parameters, 
those of the corresponding physical ones in the suitable renormahzation scheme that has 
been adopted (in our case, we are using systematically the minimal reduction scheme). 



III. RESUMMATION OF SUBLEADING LOGARITHMS IN GENERAL SUSY- 
PROCESSES 

The size of the one loop SUSY Sudakov corrections, investigated in Refs. 
Q], indicates that at TeV energies one must in principle also include the higher order 
contributions as mentioned in the introduction. While for S M p rocesses a lot of attention has 
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20 



21 



22 



25 



26 



been devoted to this particular problem in Refs. 0, [4 , 
for SUSY Sudakov logarithms only higher order corrections to scalar production in e~^e 
collisions are known Il3ll. 



We therefore discuss the situation in the MSSM for final fermion production in some 



more detail in order to clarify the arguments. As in Ref. 



we assume a "light mass" 



SUSY scenario (with the SUSY mass scale rris ^ 
for all particles involved in the loop corrections. 



rrit ~ rriH = M to logarithmic accuracy) 



For what concerns the DL and angular dependent logarithmic corrections, since they 

treatment will be identical 
so that their exponentiation 



are only mediated by the exchange of SM gauge bosons, their treatment will be identical 



with that already given in the Standard Model case 
will be granted at the same subleading logarithmic order accuracy, and does not need to be 
rediscussed here. 



For the subleading Sudakov corrections of the universal, i.e. process independent type, 
we have now novel contributions in the MSSM, both from particles with electroweak gauge 
couplings and from particles with coupling of Yukawa type. Since the former where already 
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discussed in Ref. [IS", we will only discuss here the novel Yukawa type Sudakov terms 
for third family quarks which are of particular interest since they contain a strong tan/? 
dependence. For the purposes of this paper we shall only consider the asymptotic corrections 
above the scale of electroweak symmetry breaking since the mass-gap contributions originate 
only from QED and are therefore in principle known P, ^ | . For realistic collider simulations 
they must, though, carefully be included via matching at the weak scale. 

In order to establish the "exponentiation" of the one loop Yukawa corrections, i.e. the 
fact that we only need to consider the situation depicted in Fig. (jT)), we need to show that 



the diagrams of the type shown in 



As we argued in Ref. 
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30 



ig.dZl) cancel each-other. 



31 



32| using the symmetric basis, this feature is ensured 



by the gauge invariance of the Yukawa sector. In order to demonstrate the cancellation in 
the physical basis let us first consider the case of right handed external top quarks. 

In the following we only need to consider soft gauge boson insertions (with loop mo- 
mentum /) since we want to generate three large logarithms at the two loop level. This 
means that the /-dependence of the loops with Yukawa couplings can be neglected to SL 
accuracy since the on-shell self energy and vertex diagrams don't produce an infrared-type 
contribution in the limit I ^ 0, k —>■ 0,pf = mf. Note that to logarithmic accuracy we can 
set pf = since all irti < rris ~ M and since, for now, we consider only the case of a 
heavy photon (A = M). The QED type corrections are included via matching as 
indicated above. We thus need to show that the UV-logarithms from the sum of the self 
energy diagrams cancel precisely the UV-logarithms in the sum of the vertex corrections. 

The generic one-loop diagrams (modulo the couplings and with common mass scale ~ 
rrit ~ M), corresponding to the inner loop insertions depicted in Fig. Q, are given by: 

c ^ r ^ 

^ ~ J (27r)"(A;2-m2)((A;-pi)2-m2) ^ ' 



(27r)" (P - ml){{k + pi)2 - m2) 

for the self energy insertion and for the vertex diagrams we have for zero momentum transfer 
("soft gauge bosons", i.e. / = 0): 

J (27r)"(A;2 -m2)((A;-pi)2 -m2)2 ^ ' ' 

2 ^ r rk + + 

J (27r)"(A;2-m2)((A;+pi)2-m2)2 ^ ' ^ 
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Now it is straightforward to obtain from Eq. ()3.ip : 

-Si (3.5) 



and analogously from Eq. ()H.2j) : 



Si (3.6) 



dpi 



/J 



Thus, for identical Yukawa couplings, we have to show that the sum of the couplings to the 
various contributing diagrams, multiplying the same UV-divergence, cancel each-other! 

Indeed this is what happens when one makes on the one hand the sum of the two loop di- 
agrams with the vertex contributions for the Yukawa terms (diagrams on the left of Fig. Q) 
and one the other hand the sum of the two loop diagrams with the self energy insertions 
(diagrams on the right of Fig. Q). This is easy to check in the case of external right handed 
top quarks with the exchange of a photon and a Z-boson. For each type of (scalar, fermion) 
virtual contribution like (G^, 6), (i/^,6), (6, x^), (G°, t), ... (t, X°) one sees that the same 
contribution e'^Qf/c^ due to the sum of the photon and of the Z-boson exchanges factorizes 
the Yukawa couplings in the case of the self energy insertions and in the sum of the two types 
of vertex contributions. In the left handed top quark case, an analogous result is obtained 
by considering all the left handed diagrams of the process. 

As mentioned above, this cancellation is a consequence of the fact that in spontaneously 
broken gauge theories also the Yukawa sector is gauged and that softly broken supersymme- 
try preserves the gauge structure. We can therefore employ the non-Abelian generalization 
of the Gribov theorem [3| (in the context of the infrared evolution equation method 0,Q|) 
as indicated in Fig. (Q). As a result one obtains as a solution of the evolution equation, for 
all orders, to SL accuracy, the exponential of the universal (DL and SL) and of the non uni- 
versal (angular dependent) SL one loop contributions listed in Section II. In addition, there 
appears a contribution arising from the implementation of the RG effect in the couplings 
of the exchanged gauge bosons. This is explicitely discused in ref.jlll] and the result was 
already reproduced in eq.(3.6) and (3.9) of ref. (IS]. Using the running expressions, eq. ()2.4j] . 
of the gauge couplings and expanding to subleading accuracy the rather involved combina- 
tions appearing in the equations of ref.|ll|, one obtains the additonal log^ term, with 
the coefficients b, in the following equation. 
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+2(6. + b,){2 log - log^ ^) + 26j"'=(log + 26r/(log ^) } (3.7) 

64^^ + ^56^ "-"'^ ^^-^^ 

The coefficients b are extracted from the coefficients c™, c^°, c'^ukawa^ ^ang Sect. II: 

V/a(/a + l)+/^) a = a,P (3.9) 



and 



b' 



ang 

«/3 15^2 



(3.11) 



where Jq denotes the total weak isospin of the particle with chirality a, its weak hy- 
percharge, mj = and m^/ =. It should be noted that the one-loop RG corrections do not 
exponentiate and are omitted in the above expressions. They are, however, completely de- 
termined by the renormalization group in softly broken supersymmetric theories such as the 
MSSM and sub-subleading at the higher than one loop order. The one loop value is the 
already mentioned RG term Eq. ()2.5p . 

In Eq. ()3.7p we use the SUSY-mass scale in the Yukawa and SL-RG terms, but not 
in the gauge terms. Using the wrong scale in the DL-type corrections would unavoidably 
lead to wrong SL contributions. 

The already existing explicit results for sfermion production can be obtained in the 
straightforward way by using the corresponding expressions of the various coefficients b. 

In order to generalize the above results to arbitrary "light" SUSY processes it is con- 
venient to work in the symmetric basis, i.e. in terms of the symmetry eigenstates. This 
is particularly important for the chargino production that we discussed at one loop in the 
previous section. In the general case let us denote physical particles (fields) by / and parti- 
cles (fields) of the unbroken theory by u. Let the connection between them be denoted by 
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/ = J2u C-^^u, where the sum is performed over appropriate particles (fields) of the unbro- 
ken theory. Note that, in general, physical particles, having definite masses, don't belong 
to irreducible representations of the symmetry group of the unbroken theory (for example, 
the photon and Z bosons have no definite isospin). On the other hand, particles of the un- 
broken theory, belonging to irreducible representations of the gauge group, have no definite 
masses. Then for the amplitude Ai-^^''"^"-{{pk}, {nT-i}', M, A) with n physical particles fi with 
momenta pi and infrared cut-off A, the general case for virtual corrections is given by 

n 

Ul,...U„ j = l 

In the following we give only the corrections for a light SUSY mass scale ~ M and for 
a heavy photon (A = M) with all \2pipk\ ^ ml, M^. In this case, we can easily work in the 
symmetric basis and give the results for these amplitudes. As discussed above and described 



in detail in Refs. 



0, 11, 13, Ell, the soft virtual and real QED corrections must be added 
by matching at the weak scale M. It should be mentioned, however, that the Yukawa terms 
are independent of the matching terms and that the Ward identities of the type ()3.5p and 
(13. 6|) now apply to the amplitudes A^"i' - ""({pfc}; M). 

Under these assumptions, we have for general on-shell matrix elements with n-arbitrary 
external lines the following resummed SL corrections in the symmetric basis: 

A^sL"-'"" {{Pk}; m.; M) =exp{± -i (^.(log^ 1^ - 2 log ^)) . . 

1 5 1 — t5 S 

log ^ + &r log ^}aib:::"'"^'"""^' """"(K}) (3.13) 

with bk,bk, bl""^ defined in eq. (|3.8ll3.1(ljl and 



C^ = ^E E ^.^^log^ (3.15) 

°" l<kVa=B,Wi ^ 

The fields u have a well defined isospin, but for angular dependent terms involving CKM 
mixing effects, one has to include the extended isospin mixing appropriately in the corre- 
sponding couplings I^'^i^ of the symmetric basis. If some of the sparticles should be heavy, 
additional corrections of the form log^ ^ etc. would be important. 
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The result in Eq. fl3.13p . is valid for arbitrary softly broken supersymmetric extensions 
of the SM with the appropriate changes in the /3-functions. Taking the SUSY-QCD limit 

^ a,,g' ^ + 1) = + l) ^ C A = 3 , I , {l g + l) = 1^(1^ + l) ^ C F = 

4/3, Uh = 0, M = \g = rrigjCj^^ = 0) of the various terms, Eq. ()3.13|1 is also valid for the 
virtual SUSY-QCD results. It should be emphasized, however, that in this case the virtual 
corrections are not physical in the sense that the gluon mass is zero and thus we would 
need to add the virtual matching and real contributions before we could make predictions 
for collider experiments, while in the SM soft QED energy cuts can define an observable 
and the heavy gauge boson masses are physical. In any case, the form of the operator 
exponentiation in color space agrees with the dimensionally regularized terms in Ref. 
for non-SUSY QCD. 

The general result in Eq. ()3.13|) agrees on the cross section level with the specific cross 
section expressions for fermion in Eqs. ()3.7p and similar ones for sfermions and charged 
Higgses, as well as with the one loop results for chargino production presented in the previous 
section. The appropriate mixing matrices C^"^ and RG-terms must be included for this 
comparison at one loop. They are sub-subleading at higher orders. 

We have now at our disposal the sub-leading expressions for all the considered final states, 
both at the one-loop level and completely resummed. Our next goal is that of comparing the 
two approximations at variable energy and verify whether and where they can be considered 
as "essentially" (i.e. at the expected one percent experimental accuracy level) equivalent 
or, in the same spirit, "drastically" different. This will be done in the forthcoming Section 4. 



IV. COMPARISON OF THE TWO APPROXIMATIONS AT SUBLEADING 
ORDER ACCURACY 

In this Section we evaluate numerically the basic observables (cross section, forward- 
backward and left-right asymmetries), whose expression is given in Appendix A, when the 
final state is ti, bb and when it is a pair of charged "genuinely" supersymmetric partners, 
i.e., sfermions, charged Higgs and charginos. For simplicity, we do not report the analysis 
for production of leptons or light quarks. In these cases, the previous analytical expressions 
apply, with the simplification that the Yukawa terms can be neglected. The numerical 
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results are comparable with the cases considered here. We compute the full effect in the 
observables in two approximate asymptotic approaches. First, we consider the complete 
set of one loop Sudakov contributions. As we have explained, these are terms growing 
like \og^ or log-^^^. Then, we compare this result with the one that is obtained by 
resumming to all orders at sublcading logarithmic accuracy. We remind that this corresponds 
to using an expression that predicts rigorously all the terms of the form log^^ {s / My^) 
and a^log^^-\s/M^). 

The deviation on the total cross section, the forward-backward asymmetry and the left- 
right asymmetry are defined as 

^<^tot = ^ - 1 (4.1) 

Ctot 

_r ^Born 

C^tot CTtot 
_ ^Born 

A^LR = ^ - %7 (4.3) 

(J tot CTtot 

where (Jtot, ^FB, ^LR arc the radiatively corrected observables. If we consider for instance 
the total cross section, in the perturbative one-loop scheme we have 

^ ^Born I ^onc loop i\ 

t^tot = CTtot + C^tot , (4.4) 

and in the resummation scheme we must use the expressions that we have described in the 

previous Section. 

For each final state and for both the one-loop scheme and the resummed scheme we 
consider two values of the important parameter tan /? that we choose to be tan (3 — 10 and 
tan P — 40. This will allow a discussion of the role of the phenomenologically important 
Yukawa terms and of the validity of the two approaches. 

The one-loop Sudakov terms are quadratic and linear (and the resummed expressions are 
based on them). At sublcading accuracy we need not specify the scale of the linear logs. The 
scale of the quadratic ones is conversely important at this level of accuracy. In conclusion, 
the choice of scales can be numerically relevant and we explain how we fixed it. 

quadratic logarithms: these are terms of purely gauge origin and their scale is pre- 
dicted unambiguously by expanding the full one loop calculation. It is Mw or Mz 
or the infrared photon mass regulator My according to which of the various diagrams 
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originating them is considered. For the aim of our discussion (mainly the comparison 
with the resummation approach) we have used Mw (called My in the Appendices) in 
all such terms. 



linear logs of gauge origin: these are the logarithms that combine with the quadratic 
ones in order to reconstruct the 2 log — log^ combination. Here, we take My = M^/ 
for the same previous reasons. 

non universal linear logs: these are single logarithms multiplying non-trivial functions 
of the scattering angle. Again, these are originated by diagrams where the correct 
scale can be taken to be My = M^. 

linear RG logs: here My = Mw- 

SUSY logs: these are single logarithms of Yukawa type originated by diagrams with 
exchange of SUSY partners. We use a common scale Msusy fixed at Msusy = 300 GeV. 
Here, the choice of Msusy is rather arbitrary, the difference being a sub-subleading 
term constant with respect to the energy. Our choice is motivated by the recent 



investigation j36| where it is shown that this value is a typical one for scenarios where 
the Sudakov expansion can represent accurately the MSSM effects at energies above 
1 TeV. 

After these preliminary remarks, we turn to the discussion of the numerical results and in 
particular, of their reliability. With this purpose, we try to sketch a preliminary qualitative 
evaluation of the expectable accuracy of the resummed expressions. For simplicity, let us 
consider just the gauge logarithms, that is the combination 21og — log^. In the one loop 
approximation, the correction is of the form: 

s s 

relative effect at one loop = 1 + 5one loop, ^onc loop = Cil ■ a(2 log — log^ T72~)' i^-^) 

where cil is a numerical constant. Exponentiating this term would produce, at two loops, 
the extra correction \5'l^^ i^^p. Since the resummation is able to identify correctly only the 
leading and subleading logs at all orders, but not the sub-subleading terms, we would have 
for the complete two-loop correction an expression of the kind: 

1 s 
real relative resummed effect = 1 + 5one loop + -^ono loop ± ^^a^ log^ j-j^ H (4.6) 
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where we have used in the theoretical error the smallest scale that we have. The dominant 
theoretical uncertainty of the resummation procedure can then be estimated at 3 TeV as 

±^2^^neloop (2 log - log^ ^ )2 " ^ 28 C?l ""^^ '""^ ^ ^ 

Hence, the dominant theoretical uncertainty of the resummation expressions (at subleading 
accuracy) would be below 1% if the one loop effect were below 50%, 25% or 10 % for 
cr/cil = 1.1,2.1,5.3 respectively. Since potentially large contributions to cr can in 
principle appear due to terms that are already present at one loop, (e.g. relatively large 
non logarithmic - constant - terms like those found in for charged Higgs production) 
and since the above analysis is admittedly naive, our conservative attitude would be for the 
moment that of not considering the resummation expressions as the final word for a high 
accuracy prediction when the one loop effect is beyond, say, 10%, in which case we feel that 
a (tough!) partial two- loop calculation would be highly desirable. 



In the remaining part of this Section, for each final states, we shall plot the effects in 
the two asymptotic approximate approaches as function of the energy for two values of 
tan/5 = 10,40. We also collect in Tab. (jH) the results for tan/5 = 10 at 1 and 3 TeV. The 
list of the considered case is now following. 

final top and bottom 

The effects are shown in Figs. (I3p4|l . At 1 TeV, the one- loop effects in at and are still, 
essentially, within the assumed limit of 10%, and their differences with respect to the 
resummed expansions are of approximately one percent. Similar conclusions apply for the 
two considered asymmetries. To the subleading logarithmic accuracy, a one-loop calculation 
seems therefore sufficient for all the considered observables at that energy, i.e. it does not 
"practically" change after (the corresponding) resummation. At 3 TeV, the one-loop effect 
on at is definitely beyond the ten percent value, the difference with respect to the resummed 
expressions reaches a ten percent size (note the reduction of the size of the resummed 
effect). Also the two asymmetries show visible differences between the two approaches (and 
^LR is particularly large at one-loop) at that energy, as one sees from the Figure. The 
situation appears definitely better for bottom production, where both the one-loop effects 
and the differences with respect to the resummed expressions remain reasonably small 
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(also, it is conceivable that in this case the available experimental accuracy is worse than 
in the top case). 

final sfermions and charged Higgs 

The effects are shown in Figs. (jSlini)- At 1 TeV all the observables are "under control" with 
quite small correction from the resummation procedure. At 3 TeV, there are problematic 
large corrections to cr^-^, Ap^^^, a^^, ai^, Ap^^^, an, Aps^n- In particular, for the cross 
sections a^^, a^^, cr^^, and an the effects are around 20 % and raise, in our opinion, serious 
doubts about the reliability of the resummation at this level of accuracy, 
final charginos 

If the final state is a pair of chargino and antichargino, we can in principle consider three 
separate cases according to which charginos are produced. In other words we consider the 
process 

e+e" ^ xtXb (4.8) 

in the three cases {a,b) = (1, 1), (1, 2), (2, 2). In this case the observables depend on the 
mixing between the Higgsino and the Wino component in the mass eigenstates charginos. 
The chargino mixing matrix depend on the MSSM standard parameters M2, /i and tan/3. 
We have chosen to evaluate the observables at M2 = 200 GeV and fi = 300 GeV as a 
representative light scenario. In a forthcoming paper, we shall discuss in full details the 
dependence on the mixing focusing on the constraints that virtual correction impose on the 
determination of the chargino system parameters. The effects are shown in Figs. ()1()II12|1 . In 
the phenomenological analysis, it is sometimes reasonable to assume that only the lightest 
chargino has been produced asking for what information can be gained in this case. However, 
especially in a light SUSY scenario, one can also consider a more favorable situation where 
both charginos can be produced. Then, one can study the inclusive observables like, for 
instance, 

O-inc = 0-11 + 20-12 + 0-22 (4.9) 

for which the Sudakov expansion is completely independent on the mixing and there is no 
dependence on both M2 and fi. The effects in this representative example are shown in 
Fig. (O. Looking at Tab. (P) we see that at 1 TeV the asymmetries have corrections below 
10% and the cross section is just above that value. The effect of resummation is in the range 
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of 1-2 percents for a and ApB and one order of magnitude smaller in A^r. At 3 TeV, on 
the contrary, both the cross section and the forward-backward asymmetry are large and the 
resummation introduces large shifts, particularly in the cross-section (-1-10%). A detailed 

analysis reveals that the large effects are mainly due to the Wino component where, at one 
loop, there is not single logarithm compensating the leading squared logarithm. 

V. CONCLUSIONS 

In this paper we have considered all the processes of production of pairs of charged par- 
ticles with spin zero and one-half from electron-positron annihilation at TeV energies in the 
framework of a given (MSSM) supersymmetric model. Assuming a (relatively) light SUSY 
scenario, with typical values for the relevant masses below a few hundred GeV, we have 
concentrated our attention on the asymptotic Sudakov logarithmic expansion of the elec- 
troweak component of the invariant scattering amplitude. We have computed at one-loop 
all the leading Double logarithms and the subleading Linear logarithms of the expansions, 
also including the linear logarithms of RG origin. To subleading logarithmic accuracy, we 
have also computed, for the same processes, all the resummed exponentiated expressions. 
To our knowledge, this is the first complete calculation of this kind, and we do not have at 
disposal calculations of different authors with which to compare our expressions. We have, 
though, performed an internal self-consistency check of our results for the chargino case, and 
we hope that there are no mistakes in our several formulae. 

Given the two approximate asymptotic expressions, we have verified that, for energies in the 
one TeV range (final goal of a future LC), there are practically never visible (in our work- 
ing assumptions, beyond a one percent level) differences between the two results for a set 
of realistic experimental observables. On the contrary, rather strong discrepancies appear 
almost systematically as soon as one approaches the few (2,3) TeV regime (goal of a future 
CLIC). 

A first conclusion is therefore that, in the one TeV range, an asymptotic Sudakov expansion 
at one loop has a subleading logarithmic accuracy that does not require extra resummation, 
for the considered processes. The same conclusion cannot evidently be drawn in the higher 
energy considered regime. 

We stress at this point two important facts. The first one is that the validity to sublead- 
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ing accuracy of the one-loop expansion does not necessarily mean that it is completely 
accurate when one moves beyond that level of accuracy. Extra terms, in particular energy- 
independent ones, might well be relevant at the one percent level in the expansion. The 
second one is that the validity of the expansion itself, i.e. the fact that it gives indeed an 
adequate description of the real (complete) effect even after inclusion of possible extra e.g. 
constant terms, remains to be demonstrated. 

It is rather simple to realize that the two points are related, and one possible solution would 
be represented by the preparation of a full one-loop program. This would allow to check the 
validity of a logarithmic expansion, at the same time allowing to determine by a proper fit 
to the complete result the value of a possible extra constant term. 

In the s pec ial case of charged Higgs production, this task has been actually carried 
through [3l|. One can see therefore that for that process an asymptotic Sudakov expansion 
at one loop, with an additional constant term to be realistically estimated, is suitable at 
one TeV. The benefit of this conclusion is represented by the fact that, in such an improved 
expansion, the coefficients of the different asymptotic terms depend each one on special 
reduced subsets of the supersymmetric parameters. This would allow, via identification of 
the various terms, stringent and simple tests of the model, of which one must assume the 
previous experimental confirmation. In fact, a declared goal of future colliders will be also 
that of performing precision tests of the (assumedly discovered ) supersymmetric model. 
In the case of charged Higgs production, the coefficient of the linear Yukawa logarithm was 
only dependent on tan [3. For the latter, particularly if its value were large (beyond, say, 
ten), the precise experimental determination is not at the moment completely clean, and we 
proposed in previous papers to measure tan /? from the slope of the final pair 

cross section. In the case of e.g. chargino production, other SUSY parameters e.g. of mixing 
type would enter the coefficients of the subleading logarithms, and different ones will appear 
in next-to subleading terms, so that another independent relevant test of the model would 
be available. To perform the test in a rigorous way would require the preparation of a com- 
plete one-loop program, analogous to the one completed for Higgs production. From what 
shown in this paper, we would hope that, in a relatively light SUSY scenario, a logarithmic 
one-loop expansion (with a possible addition of an extra constant term) were able to provide 
an adequate description of chargino production in the one TeV range, and no need of hard 
two-loop calculations were advocated. We are already working along that direction. 
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APPENDIX A: OBSERVABLES FOR PRODUCTION OF FERMIONIC AND OF 
SCALAR PAIRS 



1. Fermionic pairs 

For any type of fermionic pair (leptons, quarks, charginos, neutralinos) we can write the 
complete amplitude as: 



ab 



with a,b = L or R and 



(Al) 



{j^Par = v{e+) u{e-) {^,Pbf = u{i) j^Pa v{j) (A2) 



The unpolarized angular distribution is given by 



da 



d cos 327rs 
and the Left-Right polarization asymmetry is obtained as 



(A3) 



da^^ , da 



d cos d cos 'd 



with 



da^^ 



dcos^ 327TS 
The Forward-Backward asymmetry is 



(A4) 



(A5) 



A 



FB 



(If /b) dcosi? 

Uf + Ib) 



dcosi? 



where = ^^dcosd and = ^\dcosd. 



(A6) 



2. Scalar pairs 



Following the notations of ref. [13[ we write the amplitude as: 

2e^ _ 

A = — t;(e+) {aiPi + afl^/j) u{e 
s 



(AT) 
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The unpolarized cross section is 



= nJ-^ (|a^p + \an?) (A8) 

ct cos w 8s 

whereas the Left-Right polarization asymmetry is obtained from eq. ()A4|) with 

'^^^'^ nJ-^ smH{\aL\'-\aa?) (A9) 



d cos d 8s 
and the Forward-Backward asymmetry is still given by dMl). 
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APPENDIX B: AMPLITUDES AT ONE LOOP 



1. Leptons or quarks 

The Born amplitudes for a final fermion pair //, with the notations of Appendix A are 

^ir = j^j-Ai^^l - mi}) - AiT = - (bi) 

Air = 4- l^f - Qf] ^IT = -4- IQf] (B2) 
cis cis 

The one loop terms, factorizing these Born terms are given in Section II. 



2. Sleptons, squarks or charged Higgs bosons 



With the notations of Appendix A and ref . |13| , the Born amplitudes for a final sfermion 
pair // are 

Horn n , il'f-slQf)9eL sjQ j + (1 - 24)J| 

= -Qf + ^-^2 = ^T^2 (^3) 

Q,,(Jhj^.'l^ ,B4) 



a 



Born _ 

R - ^/ "1" 9 2^2 ~ ^2 



This applies to H^H using Qf = —1, /| = — | and Yf = —1, i.e. 

The one loop terms, factorizing these Born terms are given in Section II. 



3. Charginos xtxj 

At one loop, following the notations of Appendix A, the amplitudes A'^j, where ab refer 
to LL, LR, RL, RR, originate from s, u, t channel contributions: 



A^- = -S?f + -Ut^ + -T^ (B6) 



^2 pi pi 

^ cab ^ r/""^ _L 
— Oj • -I • -t- — 

S U t 

with 

riab Qofe, Born . qab, Born in , Qufe, Born fin . qab, ang , riafc, RG (Ti7\ 

'^ij — '^ij '^ij + Z^^ik '^kj + ~^ '^ij ' ) 
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k 

J.LR ^ J.LR, ang ^gg^ 

in which we have specified the contributions of the Born terms (photon and Z exchange in 
the s channel, sneutrino exchange in the u channel), the universal corrections from initial 
and final lines, the angular dependent corrections and the RG corrections. To make them 
explicit, it is convenient to separate the Higgsino and the Wino parts: 

At Born level, the only non vanishing terms are: 

qLL ,Hig Born 1 '7+* '7+ qLL ,Win Born 1 y+* '7+ /tdi 1 \ 



jjLL ,Win Born ^ /7+ /'T:>10^ 



1 1 

tLR ,Hig Born rv— ry—* qLR ,Win Born 



QLin, ,nig rsorn y— rv—* qLiH ,win tsorn ^ rv— ry—* ITI^ Q^ 



gRL ,Hig Born ^ _ ^ (BM) 
^Jii ,Hig Born ^ _ -Lz'^Z^-* (B15) 

such that, at one loop, one can write: 

^a.,Hig ^Igab ,Hig Born^ ^ ^ + J^f^n,H^g^ (3 log ^ - log^ ^) 

+6^ (log ^) + 6r'^^^(l0g ^) } + ^5!y-Hig,RG (B^g) 

^a6,Win ^ [^5j.WinBorn^^^.,WinBorn]^^^^^.„^^21og^-log2^) 

+i,^n,mn(_ i^g2 ^ i,-P.H^-(log } (B17) 

with the initial and final "gauge" coefficients 

Ibnsid And 
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b^'""' ^'"^ = (B19) 



the final Yukawa corrections 



WW w w 



the angular dependent terms 



C/' ^^'^ = ^r^f' = - ^ { log(7) + (1 - 7) l°g(— ) } (B23) 
and the RG terms (from the Higgsino components) 

« \7+*7+M( 1 w 3-64 + 8< 



-^w^w 



(B24) 



S - - 4^[^2.^2, ](^)( ^ ATf,^ 4 )Hl°g 



W W 



(B25) 

= - ^[^2r^5](^)(^Ar,^ + ^)[log ^] (B26) 

w 

= - £[^2.^.7](^)(^^fa. + i)[log ^] (B27) 

It is interesting to check the expression of the "gauge" Wino contribution eq. ()B19|) . by 
making the sum of the sphtting function term and of the Wino Parameter Renormahza- 



tion term, ref. 



2c£^ = [ ^] [(iVfan. - I)] [log ^] [Zt:Zt,PL + Z^.Z^ Pn] (B29) 

w 

and observing that all single log cancel and that the total Wino part is a pure quadratic log, 
exactly like in W^W^ production. 
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APPENDIX C: RESUMMED AMPLITUDES 



The case of ordinary fermions has already been given in Section III, eq. ()3.7|) . The exten- 
sion to sfermions and Higgses is straightforward using the corresponding expressions for the 
various coefficients b. We make now exphcit the case of Charginos which is more dehcate. 

In correspondence with eq. ()B10|) . following eq. ()3.13p . the resummed amplitude called B"'^, 
with ab = LL, LR, RL, RR, is written: 



= Bt-'""'^ + Bt-'^'"" (CI) 



Bo,^ = ^^St^ '"'^ exp { [U: + 6/-' ^^^] [i log=^(^)] 



+ + &^^'^'^^^)(2 log — - \oi j^) + C^log ^) + bZ'^'^'^log j^) } 



+Blf (C2) 



2 2 1 

rjab.Win _ nab ,Win Born , ^ jjab ,Win Bornn r rTm , I/in, VKm] r l^„3/^_\i 

(C)(21og^-log^^) + 6/^'^'^^'^(-log^ ' ' 



^^afe, W^n, RG ^^^^ 



with the new quantities (not defined in the one loop expression): 



Tfin, Win _ Q^^/^O j^Win PR _ '^(^0 /q^\ 



w 



4:ssi(s)ci(s) Assid ■' 



eHs) , 



BLR, m,, RG ^ _^ yy _ }lz,,Z^*] (C7) 

4ssi(s)c4(s) 4:sslcl ■> 
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s-'""-°--{^-Is^1bo...}k-^;i (C8) 



git, ,«„, ^ _j ^!(f) ^ ^!(£) _ 1^ + }\Zt-Zt,] (CIO) 



j^ilL, Win, RG _ j^RL, Win, RG _ q (C12) 

where the running expressions depending of s are defined in eq. ()2.4p . 

One can check that with the expansion exp(X) = 1 + X one recovers the one loop contri- 
bution apart from a two loop log^ term and a reshuffling of the Wino PR terms. 
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FIG. 1: Higher order corrections to the original on-shell vertex diagrams with Yukawa couphngs. 
The non-Abehan generahzation of the Gribov theorem can be apphed as is shown in the text. The 
figure is only schematic since in principle the gauge bosons couple to all external legs. 
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FIG. 2: Two loop Feynman diagrams involving Yukawa couplings. In the text it is shown that the 
sum of all such contributions with self energy subloops are canceled by the sum of the corresponding 
vertex diagrams to SL accuracy. 
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TABLE I: Summary table for the effect at tan (5 = 10. The numbers in boldface are effects larger 
than about 10%. 
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